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ABSTRACT 


The concept of nonlinear coherence employed in car- 
rier tracking to improve telecommunications effi- 
ciency is disclosed. A generic tracking loop for a co- 
herent receiver is shown having seven principle feed- 
back signals which may be selectively added and ap- 
plied to a voltage controlled oscillator to produce a 
reference signal that is phase coherent with a received 
carrier. An eighth feedback signal whose nonrandom 
components are coherent with the phase detected and 
filtered carrier may also be added to exploit the side- 
band power of the received signal. A ninth feedback 
signal whose nonrandom components are also coher- 
ent with the quadrature phase detected and filtered 
carrier could be additionally or alternatively included 
in the composite feedback signal to the voltage con- 


trolled oscillator. 
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4 Claims, 6 Drawing Figures 
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COHERENT RECEIVER EMPLOYING 

NONLINEAR COHERENCE DETECTION FOR 
CARRIER TRACKING 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958 Public Law 85-568 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to coherent receivers and em- 
ploys the concept of the nonlinear coherence of ran- 
dom nonlinear oscillations, and move particularly to re- 
ceivers in which nonlinear coherence is employed to 
increase telecommunication efficiency. 

The meaning of the term “nonlinear coherence” can 
be explained as follows. The usual method for examin- 
ing the mutual power between two signals s,(t) and 
j 2 (t) at an arbitrary frequency /is through the use of 
the so-called cross-spectrum. In essence, the cross- 
spectrum of two signals will represent the spectral den- 
sity of power that is mutually shared in a phase coher- 
ent manner. It is important to note that each signal can 
have power in the same frequency band without there 
being cross-spectral power in that band. Thus, having 
common frequency components does not guarantee 
mutually coherent power. On the other, in order to 
have mutually coherent power, the signals must be 
phase coherent. 

In a nonlinear system, it is possible to have coherency 
between a signal at one frequency, say f l} and another 
signal at some multiple of that frequency. For example, 
in a squaring loop which is commonly used for tracking 
suppressed carrier signals, such an input signal with en- 
ergy centered around f\ is squared (nonlinear opera- 
tion ) to produce a signal centered around 2 f x which is 
phase coherent with the suppressed carrier signal. The 
signal at 2/, is then tracked by a conventional phase- 
locked loop with a VCO whose nominal frequency is 
2 /,. 

Mathematically speaking, consider the nonlinear 
cross-spectrum 

Su(/; m, n)=J'J_^s 1 “(<)*j n (<+r) exp (— i'2ir/r)#dT 

( 1 ) 


where m and n are integers. To test for the possibility 
that s 2 (t) has a frequency component coherent, in a 
nonlinear manner, with a component in Ji(t) at twice 
the frequency, one would compute 

Su(/; 2, l)=JJ_ aj 8i a (t)s J (t+T) exp ( — 2iri/T)dtdT 

( 2 ) 


In terms of the squaring loop, j ( (t) would be the bi- 
phase modulated (suppressed) carrier input signal and 
j 2 (t) the phase-locked loop reference signal at 2/j, i.e., 

s,(t) = d(t) sin (2nf x t) 

j 2 (t) = cos [2rr(2/i)t] 3. 


2 

If for the moment we ignore the modulation d( t) on 
s,(t), then Si 2 (f;2,l) would have a spectral line at 2f x . 
Thus, the signals r,(t) and r 2 (t) are said to be nonlin- 
early coherent and the coherent receiver structure of 
5 the present invention is conjectured on this principle. 

As satellite and deep space technology have ad- 
vanced rapidly, even in the few years of its history, top- 
ics of increasing current interest are the application of 
10 Earth Satellites to the development of tracking and 
data-relay satellite networks for relaying earth resource 
data, earth-orbiting manned space/base stations, tacti- 
cal communications satellite systems, integrated 
communications-navigation networks, air traffic con- 
1 5 trol systems, etc. Outside the application of satellites in 
orbit about the Earth, interest centers around the plac- 
ing of communication satellites in orbit about Mars, 
and the sending of exploratory spacecraft to Jupiter, 
Neptune, Saturn and Pluto. While such applications 
20 impose autonomous operation of long periods of ser- 
vice on both man and machine, they also place in- 
creased demands on telecommunication system effi- 
ciency. Telecommunication system efficiency means 
the effectiveness with which a system performs both 
25 the tracking and the communication functions. In what 
follows we develop the theory as it applies to the vari- 
ous areas of carrier and suppressed carrier tracking, 
subcarrier tracking and phase-coherent communica- 
tions. 

30 

SUMMARY OF THE INVENTION 

In a receiver chann el fo r a time varying signal x char- 
acterized by r = V 2 P c sin 4> + V2S X cosd> + n, 
where x = 3d is a biphase modulated subcarrier, Sand 
35 d represent the data subcarrier and the data waveforms, 
respectively, which are assumed to be square wave- 
forms, and where <I> = 0 characterizes modula- 

tion due to receiver motion or the randomness of the 
channel, P c =m 2 P represents power at the carrier fre- 
40 quency, S = ( 1— m 2 )P represents the power remaining 
in the modulation sidebands and m denotes the modu- 
lation, and where Sjand A represent the receiver’s esti- 
mates of the data subcarrier and the data waveforms, 
respectively, a generic tracking loop, provided to ex- 
5 ploit the principle of nonlinear coherence is comprised 
of: a voltage controlled oscillator for generating a time 
varying reference signal r u = VT K , cos $; a summing 
junction and a smoothing filter coupling the junction to 
a control terminal of the oscillator; a 90° phase-shift 
u network for providing a quadrature phase reference 
signal r, = \T2 A 2 sin $ where $ is the time varying 
loop estimate of 4>; two multipliers responsive to the 
receiver signal and the signals r u and r, for producing 
55 quadrature phase error signals e„ —xr u and f.f=xr t \ a first 
low-pass filter of a particular bandwidth and gain cou- 
pling the signal r u to a point(2)connected to the sum- 
ming junction; a first multiplier having one terminal 
connected to receive the output of the first filter and 
60 the output of a second low-pass filter of a particular 
bandwidth and gain to provide a product signal at a 
point(l)connected to the summing junction; means for 
demodulating the phase error signal e u by a phase esti- 
mate of a reference square-wave subcarrier and a third 
65 low-pass filter of a particular bandwidth and gain for 
filtering the demodulated signal; means for delaying 
this subcarrier demodulated and filter signal a time T 
equal to a data symbol period; means for multiplying 
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this first delayed signal by a ((— T) where d(i) is the 
time varying estimate of the data waveform; means for 
multiplying the output of this last multiplying means by 
the output of the first filter to produce a third feedback 
signal at point(3)connected to the summing junction; 
means for demodulating the phase error signal e t by a 
phase quadrature estimate of the reference square- 
wave subcarrier and a fourth low-pass filter of a partic- 
ular bandwidth and gain for filtering the phase quadra- 
ture demodulated signal; means for delaying this sub- 
carrier phase quadrature demodulated and filtered sig- 
nal the period T; means for multiplying this second de- 
layed signal by d(t— T) to produce another signal at a 
point@connected to the summing junction; means for 
multiplying the output of this last multiplying means by 
the product of the multiplying means of the first de- 
layed signal and <J(r— T) to produce yet another signal 
at a point(5)connected to the summing junction; means 
for multiplying the output of the second low-pass filter 
and the output of the penultimate multiplying means to 
produce a signal at a point®connected to the summing 
means; fifth and sixth low-pass filters of particular 
bandwidth and gain connected to the outputs of respec- 
tive third and fourth filters; and a multiplier having its 
output terminal connected to a point®, one input ter- 
minal connected to the output of the sixth filter and an- 
other input terminal connected to the output of the 
fifth filter by an operator which provides a function ap- 
proximately equal to tanh x, where x is the output of 
the fifth filter. The gain of these filters may be selec- 
tively set to zero to effectively remove signals at points 
®through(7)to provide a desired combination of feed- 
back signals to the summing junction, as for an adapt- 
ive filter, or to optimize tracking for a particular appli- 
cation with minimum hardware, in which case circuitry 
associated with only disconnected feedback signals 
may be omitted. 

To exploit sideband power in applications where 
phase error can be assumed to be constant over several 
data symbol intervals, additional feedback signal at 
points(§)and(9)may be provided by a delay means of a 
period T coupling the receiver input signal x to two 
multipliers receiving the reference signals r u and r t , sep- 
arately for phase detection of the delayed input signal, 
and separate low -pass filters of particular bandwidth 
and gain coupling the outputs of the multipliers to the 
points®and(9). The filtered output of the inphase error 
signal thus produced at point®is connected to the sum- 
ming junction, and the filtered output of the quadrature 
phase error signal thus produced at point® is con- 
nected to the output of the second filter for addition to 
the corresponding quadrature phase error signal fil- 
tered through the second filter. The nonrandom com- 
ponents of these signals at points(§)and(9)are coherent 
with the corresponding signals at the outputs of the first 
and second filters, but their noise components are or- 
thogonal in time with those corresponding signals. As 
in the case of feedback signals at points(T)through®, 
the feedback signals at points® and® may be selec- 
tively removed, either actually or effectively by reduc- 
ing their filter gain to zero. However, the feedback sig- 
nal at point(8)is advantageously connected only when 
the feedback signal at point® is connected and the 
same is true for points(3)and® . 

The feedback signals at points® through® may be 
advantageously provided in all possible combinations 
taken 1, 2, 2, 4, 5, 6 and 7 at a time, and the feedback 


4 

signals at points®and®may be added to these combi- 
nations to form additional combinations with the limi- 
tations expressed or implied with respect to these last 
fe edback signals. All combinat i ons are n ew except (T), 
5 ®7 and® individually, the combination of signals at 
points® and®only, and the combination of signals at 
points® and(4)for low signal-to-noise ratio where the 
operator provides the function tanh x = x, jc«1, for 
low signal-to-noise ratios. For high signal-to-noise ra- 
10 tios, tanh x ~ sgn x, x» 1 , to provide a new combina- 
tion of just the signals at points®and(4). 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram of a generic tracking 
loop according to the present invention. 

FIG. 2 is an addendum to the circuit to be added in 
particular cases to the generic tracking loop of FIG. 1. 

FIG. 3 is a schematic diagram for an arrangement to 
0 be used to develop the signal d(t— T) in FIG. 1 for par- 
ticular cases. 

FIG. 4 illustrates a particular case of the generic 
tracking loop, namely a modified data-aided tracking 

25 lo °P' 

FIG. 5 illustrates another particular case of the ge- 
neric tracking loop, namely a modified hybrid loop. 

FIG. 6 illustrates the combined data-aided loop of 
FIG. 4 with the hybrid loop of FIG. 5. 

30 DESCRIPTION OF THE PREFERRED 

EMBODIMENTS 

A generic tracking loop for a coherent receiver 
which fully exploits the principle of non-linear coher- 
35 ence is shown in FIG. 1. The receiver is novel in that 
it suggests adding one or more levels of technology to 
that which exists in present-day tracking and communi- 
cation receivers. It also provides for the planning of fu- 
ture tracking and coherent communication systems. 
40 Many special cases of the general structure exist. 

This receiver is concerned with only a single channel 
system where the random oscillations of the received 
signal can be ch aracterized by 

x(t) = V 2 P sin [to 0 /+(cos~ , m) X(f) + 0(f)] + nt(b) 
45 

where x(t) =d(t )5(f) is a bi-phase modulated data sub- 
carrier, 0(t) characterizes the modulation due to vehi- 
cle motion or the randomness of the channel, and « ( (t) 
is a narrowband, “white” Gaussian noise process of 
50 double-sided bandwidth W ( Hz and single-sided spec- 
tral density N„ watts/Hz, i.e. 

n,(l) = VT [n<-(t) cos(a»o/+0(/)) -«,(/) 

sin(a>of+0(/))] 5. 

55 The noise components n f (t) and n,(t) are independent 
"white” Gaussian processes with single-sided spectral 
density N 0 watts/Hz. Also in Equation (4), the parame- 
ter cos“'w represents the system modulation index and 
m denotes the modulation factor. The generalization of 
u that which follows to N channel systems is straightfor- 
ward in view of that already presented in W. C. Lind- 
sey, “Determination of Modulation Indexes and Design 
of Two-Channel Coherent Communication Systems”, 
IEEE Transactions on Communication Technology, 
Vol. Com- 15, pp. 229-237, April, 1967, and W. C. 
Lindsey, “Design of BlockCoded Communications Sys- 
tems”, IEEE Transactions on Communication Tech- 
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5 6 

nology, Vol. Com-15, No. 4, pp. 524-534, August, will be described more fully hereinafter, this assumes 
1967. the phase error is constant during the symbol time T, 

The data subcarrier 5 (t) is assumed to be a square i.e., = <f>(t— T). When multiplied by the PLL feed- 

wave, i.e., a sequence of ±l’s occurring at the subcar- back signal at point(2)through a multiplier (detector) 


rier rate and the data sequence d( t) is also character- 
ized by a sequence of ±1 ’s occurring at the symbol rate. 
Extension to the case of sinusoidal subcarriers follows 
the approach taken in the last reference cited, supra. 
Assume that the ±1 ’s in the data sequence occur with 
equal probability and have a duration of T seconds. 
Under these assumptions Equation (4) can be rewritten 
as 

x = V 2 P c sin<t> + V2S X cosd> +Mj 

= a) 0 t + 6 6. 

where P c = m 2 P represents the power which remains at 
the carrier frequency and S = ( l—m 2 )P represents the 
power remaining in the modulation sidebands. If m = 
0, we have complete suppression of the carrier and if 
m = 1 we have no power in the modulation sidebands. 
In Equation (6), note that the time variable has been 
suppressed by letting x(t) = x, 6(t) = 6 , X(t) = X, <£(/) 
= <f>, etc. This will be convenient throughout for discus- 
sion, and in the drawings, although in practice it is evi- 
dent that the„signals are time varying. 

In FIG. 1, 5 and a! respectively represent the local re- 
ceiver’s estimates of the data subcarrier and the data 
waveforms. The received signal jc is applied to a multi- 
plier (detector) 10, such as a double balanced diode 
mixer, having its second input connected to a voltage 
controlled oscillator (VCO) 11. The output of the de- 
tector is connected to a low-pass filter 12 designated 
LPF, with a gain g, to indicate a low-pass filter of a par- 
ticular bandwidth and gain. The output of filter 12 is, 
or may be, connected to a summing junction 13 
through points(5)as indicated by a dotted line. Each of 
the connecting points to the summing junction repre- 
sented by a small circle at the end of a line is or may 
also be connected to another point in the circuit repre- 
sented by a small circle and having the same number in 
the circle, as shown for the connection of points (2) 
which provides conventional phase-locked loop (PLL) 
feedback to the VCO through a smoothing filter 14. To 
this basic PLL, additional elements are added as shown 
using filters of designated d-c gain (g„), where the gain 
may be zero, i.e., where the filter may be an open cir- 
cuit and the circuit between the filter of zero gain and 
the summing junction may be omitted. 

A 90° phase shifter 15 couples the output of the VCO 
to a multiplier (detector) 16 to produce a phase error 
signal €/ in phase quadrature with the phase error signal 
e, out of the detector 10. When passed through a low- 
pass filter 17 of the same particular bandwidth and gain 
as the filter 12, a feedback signal is produced which, 
when multiplied with the output of the LPF 12 in a mul- 
tiplier 18 and fed back to the VCO through the sum- 
ming junction, provides a special case of what may be 
referred to as an N-Phase Costas (I-Q) Loop, where 
N=2. See “Carrier Synchronization and Detection of 
Polyphase Signals”, IEEE Trans., Vol. Com-20, No. 3, 
June, 1972, pp. 441-454 at pages 447 and 448. 

The phase error signal e„ is demodulated by a phase 
estimate of the reference squarewave subcarrier 
through a multiplier (detector) 19 and filtered through 
a low-pass filter 20 designated LPF 2 of a particular 
bandwidth and gain, g 2 , The filtered signal is then trans- 
mitted through a T-second delay element 21 and multi- 
plied by d T = <?(/— ' T) in a multiplier (detector) 22. As 


5 23, a feedback signal to the VCO is produced at point 
d)and added to other feedback signals. 

If the output of the filter 20 is further filtered by a fil- 
ter 24 designated LPF 3 of gain g 3 and multiplied by a 
phase quadrature signal developed similarly through 
10 elements 24, 25, 26 and 27 in a multiplier (detector) 
28, a fourth feedback signal is produced at point@and 
added to other feedback signals. An operator Q( ) is 
introduced by an element which as described with ref- 
erence to Equation (7), infra., which for high signal-to- 
15 noise ratios (high SNR’s) is sgn(x). For low signal-to- 
noise ratios (low SNR’s) the output of operator Q( ) is 
simply x, where x represents the signal at the output of 
the filter 24, and not the input signal to the tracking 
loop. 

20 A fifth feedback signal at point(5)is produced by mul- 
tiplying in a multiplier (detector) 33 the output of the 
multiplier 22 by a phase quadrature signal similarly de- 
veloped through elements 31 and 32. The phase quad- 
rature signal developed in that manner at point ©is 
25 added to other feedback signals, and multiplied with 
the output of filter 17 in a multiplier (detector) 34. The 
product at point(7)is added to other feedback signals. 

In the circuit just described, the operator Q( ), linear 
30 or nonlinear, is inserted for the sake of generality. In 
practice, it is determined by the design engineer whose 
choice is influenced by the theory of continuous non- 
linear filtering. Based upon the method of estimation 
described by S. Butman and M. K. Simon, “On the Re- 
35 ceiver Structure for a Single-Channel Phase-Coherent 
Communication System,”- JPL Space Programs Sum- 
mary, Vol. Ill; No. 37-62, pp. 103-108, and J. J. Stif- 
fler, “A Comparison of Several Methods of Subcarrier 
Tracking,” JPL Space Programs Summary, Vol. IV., 
40 No. 37—37, pp. 268-275, one might set Q(x) = 
tanh(jr), although from the point of view of continuous 
nonlinear filtering theory this choice is a suboptimum 
one. Nevertheless, if such an operation were to be in- 
serted in the system, one would probably wish to imple- 
5 ment it only in one of two forms depending on the data 
signal-to-noise ratio. Since 

( sgn x x>> I 

( 7 ) 

X x<<\ 

one would remove this nonlinearity for low data signal- 
to-noise ratios and for high signal-to-noise ratios in the 
data stream, one would mechanize it by a hard limiter 
55 characteristic. The details which motivate such a non- 
linear structure will be elaborated on in what follows. 

The oscillations r u appearing at the input to the upper 
phase detector 10 are characterized by 
60 r u (t)= Vc~2 X, cos<t>(/) 8. 

while the oscillations r, appearing at the input to the 
lower phase detector 16 are characterized by 
r,(t) = VT K 2 sin<f>(t) 9. 

65 where is the loop estimate of 4>. Before proceeding 
with the derivation of the stochastic integro-differential 
equation of operation for the multiple loop configura- 
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tion of FIG. 1, one additional concept will be briefly in- 
troduced because it can easily be carried along in the 
analysis which follows. 

For a great many applications (e.g., medium to high 
rate telemetry) the phase error <t> = 6 — 9 can be as- 
sumed to be constant over several symbol intervals; 
hence, delay elements such as those in FIG. 1 can be 
used to exploit the sideband power in much the same 
manner as is done in differentially coherent detection 
or time diversity reception. An example of how this 
might be done is illustrated in FIG. 2 using elements 35 
through 39 where the phase error is assumed to be con- 
stant for T seconds and the correlation time of the addi- 
tive noise is much less than T. Elements 36 and 38 cor- 
respond to respective elements 10 and 16 of FIG. 1, but 1 5 
are in addition to and are connected to receive inde- 
pendently the reference signals r u and r,. The input sig- 
nal x is applied directly to the delay element 35 in addi- 
tion to the elements 10 and 16 of FIG. 1. In effect, two 
signals are produced at points(g)and(§) whose nonran- 20 
dom components are coherent with, but whose noise 
components are orthogonal in time with, the corre- 
sponding signal components at the outputs of the two 
filters 12 and 17 in FIG. 1. Thus, for example, one 
might add the signal at point(S)into the multiple sum- 
ming junction 13 and/or add the signal at point®to the 
output of filter 17 before further processing in the loop. 

In practice, both would normally be included together, 
or both omitted. However, the signal at point(8)is ad- 
vantageously connected to the summing junction only 30 
if the signal at point®is connected. In the ideal case, 
including them would improve the signal-to-noise ratio 
at each of these points(f)and(§)by 3 db. 

A mathematical description of the signals at points® 

— (§)will now be presented, in each case indicating which 35 
are individually similar to present day telecommuni- 
cation system designs, and which are novel. Collec- 
tively, in all possible combinations, except just the sig- 
nals at poinW®®,®» and(Jpindividually, and the com- 
bination of points® and®, and the combination of 40 
points®and®for low signal-to-noise ratios, they are all 
novel. A stochastic integro-differential Equation (26), 
infra., governs the operation of a loop which uses all of 
these signals as sources of coherent energy for im- 
provement of telecommuniation efficiency. A loop 
which uses all feedback signals might not necessarily 
yield the best performance for all applications. Only 
after a given application is analyzed will one be able to 
specify which or what combination of the signals®-® 
should be used. The generic tracking loop described 
with reference to FIGS. 1 and 2 provides for the most 
general system based upon the principle of nonlinear 
coherence. Some examples will be given which are spe- 
cial cases of the general system. However, it should be 
understood that the present invention is not limited to 
those examples. In this sense, the paper should be 
looked upon as presenting some new ideas but not an- 
swering all questions relative to their application. One 
skilled in the field of communication system theory and 
well acquainted with the published literature on the 
subject should not find difficulty in applying the ge- 
neric tracking loop to suit his particular needs by effec- 
tively selecting a gain of zero for some filters by omit- 
ting them together with signal components that follow. 
One may even find it advantageous to include all filters 
and the signal components that follow in order to pro- 
vide for switching the gain of some filters to zero under 


certain conditions, i.e., to provide for mechanization of 
an adaptive tracking loop in a coherent receiver. 

How the concept of coherence of random nonlinear 
oscillations can be exploited to the advantage of the 
telecommunication engineer by this invention will now 
be presented. We begin by presenting the equations ’ 
which represent the random voltages appearing at 
points one® through seven® in FIG. 1. 

Neglecting double frequency terms, the output of the 
upper phase-detector 10 is given by 
e u = Kj [ y/~F c sin</> + VS X cos<l> + n„(/,4>) ] 10 . 

while the output of the lower phase-detector 16 is 
= /w 2 [ v7^cos(/>— >/5 X sin<t> - n,(t,<t>) n. 

where 

n u (t,<f>) = n c (t) cos cf> — n,(t ) sin <t> 

«i(r,<|>) — n c (t) sin 4> — n,(t) cos </> 12 . 

and all multipliers are assumed to have a gain of unity. 
Clearly, n u (t,<i>) and n,(t,4>) are uncorrelated. Further- 
more, if the loop bandwidth is narrow relative to W f , 
then they can be approximated by statistically indepen- 
dent low pass “white” Gaussian noise processes of sin- 
gle-sided spectral density N 0 watts/Hz. Assuming that 
the low-pass filters 12 and 17 do not pass the modu- 
lated data subcarrier components, then the random 
nonlinear oscillations appearing at point®are given by 


o . , \KiKiP a sin 2<t> 

g “ 


+yfPllK, cos 07i„i — K t sin 0ni,] — nainnj 


(13) 


45 


50 


55 


60 


65 


where /?„, and n n are respectively the noise processes 
which emerge from the filters 12 and 17, and g, is the 
d-c gain of these filters. These processes are approxi- 
mately independent, low-pass band-limited and have 
spectra determined by the passage of “white” noise 
through the normalized filters, i.e., 

S nal (*>)=S„ (w) = /V„lG, (o>)/G,(3220)l 2 /2 

where G,(w) is the transfer function of the filters. Also 
note that g, = G,( 0). 

Neglecting double frequency terms the signal at point 
two,®, is given by 

%)= gi Ki [ vTc sin</> + n ul ] ( H) 

This signal represents the dynamic phase error in con- 
ventional PLL tracking receivers. 

Referring now to the signal which appears at point 
three®of the loop, assume for the moment that the ref- 
erence squarewave subcarrier Sis perfect, i.e., h = 5- 
This is not too restrictive since this is largely true in any 
efficient coherent receiver. Since the output of the 
upper low-pass filter 20 of d-c gain g 2 , can be repre- 
sented by g 2 K,[d VTcos <jri-n u2 ], then the delayed ver- 
sion when multiplied by d T and the signal at point®pro- 
duces 

+ VSdxdT cos -\'d T n^in^3T ] (15) 

where we have assumed that the phase error is constant 
during the symbol time, i.e., <#>(/ — ) = </>(f— 7~) and as 
previously mentioned the T subscript denotes a T sec- 
ond delay version of the corresponding signal, e g., d T 
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— d(t-T). In applications where this is not the case $§> ogy, Vol. Com-20, No. 1, February, 1972, pp. 53-55. 
does not hold since <£(/) ** <5(r— T). The spectral den- The signal appearing at point five®is given by 
sity of the low-pass approximately Gaussian noise pro- _ 

cess n u 2 is given by S-. (w) = /V o |G 2 (a>)/G 2 (0)| 2 /2 S^g^KiK^irf. d T VS cos q-f-nusx)^T 


where G 2 (oj) is the transfer function of the filter 20. 5 
Also, is approximately independent of n ul and n n 
since its energy comes from a narrow-band region of n u 
centered around the subcarrier frequency. 

When the phase-error is constant for several symbol 
intervals, tfyen d 7 d T can be replaced by its statistical av- 10 
erage E(djd r ) = 1 — 2P fi (</>) and Equation (15) reduces ■■ 
to 

s®=mA , [Vsp:(i-2p e (^)) 

+-</PldT sin 0no3T+ViS(l — 2P E (q)) cos 


(dxVSsinq + nuT ) ] 

=g 1 1 gig 1 [ ‘ S Sl ° 2 * +d T VS (sin ^ 

+ cos +nu>TnurJ 

( 20 ) 

This signal also represents energy in the sidebands 
which is coherent in a nonlinear way at the carrier fre- 
quency. 

The signal at point six(§)is given by 

%i=g2^2[ sintfrWrWfjT-] (21). 


+ rfT»ulWl!3T 


Assume that d is obtained by a matched filter technique 20 
as in FIG. 3, where an integrator 40 is followed by a 
sample and hold circuit 41 to hold the output of the in- 
tegrator at the end of an interval T until the next inter- 
val, where the interval is established by a symbol syn- 
chronizing signal of the receiver employing the present 25 
invention. A hard limiter 42 follows the sample and. 
hold circuit to yield the signal cl(t—T). The integrator 
is reset at the end of each time interval T. Then, for the 
special case of phase-shift keyed signals, the function 
Pe(4>) is the conditional symbol error probability given 3 
by ... 


where the phase-error is again assumed constant during 
a symbol time. This signal component arises in the 
data-aided loop described by the inventors in “Data- 
Aided Carrier Tracking”, IEEE Trans., Vol. Com- 19, 
No. 2, April, 1970, pp. 157-168 and in U.S. application 
Ser. No. 101,354, filed Dec. 24, 1970. 

The signal at point seven(7)is analogous to the signal 
at point<3)with the phase error 4> shifted by 90°. It is 
given by 

SQ=gig2Kf(y/P~ c cos 

(VS(1— 2Pe(*» sin -f" ) 

= gig*Kf [ VP5( 1 - 2P b U) ) 


(17) 

where R = ST/No. 

The signal appearing at point four@will be charac- 
terized for two conditions, viz., for high and for low sig- 
nal-to-noise ratios. We assume, without loss in general- 
ity, that Q(x) = tanh x. For high data stream signal-to- 
noise ratios tanh x ~.sgn x and , 

%= gigiK* [ V ? k d sin<f> + d n, 3 ] <I8>. 

i 

where d represents the data stream “estimate pro- 
duced by the upper filter 20 in cascade with sgn x , and 
g 3 is the d-c gain of the filters 24 and 27. The noise pro- 
cess is approximately low-pass Gaussian and has a 
spectral density S a ,*_=N 0 | G 2 {<d)G 3 (<i))IG 2 (0)G 3 (o) | 2 /2 
where G 3 (w) is the transfer function of the filters 24 
and 27. For low data stream signal-to-noise ratios, tanh 
,r = .r and 


+ VPo^t cos $nuT — d-rrinnm — V<S( 1 — 2Pe(</>)) sin <j>nu 


The signal at point(§)of FIG. 2 is easily found to be 
! 3D=giKi[ sin<f> + n UIT ] (23). 

40 where it is again assumed that the phase-error is con- 
stant over the symbol interval. The process n ulT is or- 
thogonal to the processes n ul , n u2 , n u3 , n tl , n t2 , n t3 , when 
the correlation time of the noise is much less than T. 
The signal at point(9)is given by 
45 S^=g t K 2 [ VT C cos«f> + n, 1T ] (24). 

when 4> is constant for T seconds. The signal S^gand S<g) 
could be added to the outputs of the upper and lower 
filters 12 and 17 of FIG. 1 to produce mutually coher- 
ent energy for tracking. 

The instantaneous phase estimate 8 or 6 which the 
receiv er pro duces is given, in operator form, by 


S^faWKiK^di/S cos 0+n u ,)(dVS sin <£+ni,)] 
=gi l gdK\K 2 £ 2 — -+dVS(cos <#>nu 4- sin <t>nai 

+*Ki« n uJ (ig) 


55 


g KyF(p) 



(25) 


where K r is the gaii) of the voltage control oscillator 
VCO. Since <t> = 8-6 we have that 


The noise process n u3 is modeled exactly the same way 
as n t3 and has the identical spectral density as n t3 but is 
approximately independent of it. Equations (18) and 
(19) represent signal energy which is mutually coher- 
ent at the carrier frequency in the sidebands and arises 65 
in a hybrid loop proposed by one of the inventors, W. 

C. Lindsey at the 1970 International Communications 
Conference in San Francisco, California, and published 
in the IEEE Transactions on Communication Technol- 


«— <»> 

This represents the stochastic integro-differential equa- 
tion of loop operation. It is the general result from 
which all loops evolve. For instance, when the symbol 
rate is such that the loop phase-error is not constant 
(typical of command systems and low-rate-coherent 
— telemetry systems) during a symbol time, the loop 
equation is given by 
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^ =g -gvm [5(D+ g (i)+ s (j)+ g (j)] (2?) 

where %>%>%» and Sj, are found from Equations (13), 
(14), (15) ana (18) respectively. 

By applying the diffusion approximation described by 
R. L. Stratonovich, Topics in the Theory of Random 
Noise, Gordon and Breach, London, England, 1967, 
the probability density function of the phase error can 
be found using the general theory given by W. C. Lind-, 
sey, “Nonlinear Analysis of Generalized Tracking Sys-j 
terns,” Proceedings of the IEEE, Vol. 57, No. 10, Octo-j 
ber, 1969, pp. 1705-172 2. Ther e it is sh own tha t 


r<t,+2r . 

p(^) = Co' exp [!/<,(<#>)] exp [— L r o(z)]dx;|<*>| <r 15 i.e.. 


Koo=J oj B[iVT(t,^)lVT(t+»’,«>)]<JT 

We also note that approximate formulas for the mo- 
5 ments of the mean time to first loss of synchronization, 
average number of slips per unit time, etc., can be 
ifound by applying the general theory given in W. C. 
[Lindsey, “Nonlinear Analysis of Generalized Tracking 
'Systems,” Proceedings of the IEEE, Vol. 57, No. 10, 
10 October, 1969, pp. 1705-1722. 

1 Some particularly interesting cases of the generic 
tracking loop will now be discussed. As noted hereinbe- 
fore, a data-aided loop is obtained by removing all 
terms from the sum of Equation (26) except Sg^nd 


Uo(<t>) = — J* ht(x)dz 


and C 0 ' is a normalization constant. For a first order 
loop H 0 (d>) is the sum of the signal terms %>through 
normalized by 2/Koo where Koo is the intensity coeffici- 
ent of the noise to be defined shortly. For low signal-to- 
noise r atios (low SNR ’s) _ 

^?Ao(<fr)=[2 ffl K,VP;+ff»«:,VS(l-2P E (0))] sin 4 

+ [gSK^S + g^KiKtS + gSK^P,] 

+ [gtg,K l ^ + g l g 1 Ki]VKS(.l-2P E {< t ,))^^ (3Q) 

while fo* high SNR the fourth term (involving g 2 z g 3 2 ) is 
replaced by the signal component g^^V - ? 
( 1-2 P t f {<}>)) sin <f> from (18). The function P e = (</>) is 
approximated by Equation (17) with R replaced by 
2S/N„W M where W 23 is the two-sided noise bandwidth 
of LPF 2 and LPF 3 filters in cascade. The above equa- 
tion represents the stiffness of the nonlinear interac- 
tions due to the multiple loops, i.e., the commonly 
called loop S-curve. 

The diffusion coefficient Kqo is determined from the 
noise and noise cross signal terms in ^through §5 For 
low SNIl’s, the equivalent total phase noise N^t,#) 
which e ffects the VC O estima te is given by 

iV T (<,0)-ffi J {VK[Kj cos 1 — Ki sin Hu] — n^nu } 

+0i'K'i(n O i+ nu,T ) 

'EfflfflKi* [fPa Sin <t>^Tfh3T 

+VS(1-2P b (4)) cos 4^01+ ^Tnuinujx] 
J rgfQi‘KiKi [dVfi ( cos 4nu+8in 4fius) 4-Wujrcu] 
T gJKiKi [(2tV*S ( titjT cos </> + Jiuir sin 4 ) 

+ rinjTn ist] + g% Ktknn 
+ffifl's^s J [VPo^T cos 0n»T— ^rniinirr 
— VS(1— 2 Pb(^>)) sin0n u ] 


4=9-— F(p)lgiKisfP' sin <t> 

P 

+ ffsKiVS(l-2P E (4>)) sin <f> + 0i K in>.i + giK^nm] 
20 ■ (33) 

The mechanization is achieved by making the gain 
in all other unused channels zero, e.g., omitting all 
other channels. In practice, the filter 12 can also be 
25 omitted in the mechanization since the loop filter 14 
will serve the same low-pass filtering purpose. 

A slight generalization of the data-aided loop is ob- 
tained by adding at the summing junction 13 the signal 
Sg)Of FIG. 2. The loop equation of operation then be- 
30 comes 

[~2 ff ,lGVPl sin <t> 

P L 

+9jjK»VS(1 — 2 P E (<f>) sin 4+9i^i(»»»i n aiT)+ffi^TnuT] 


The mechanization of this modified data-aided loop is 
illustrated in FIG. 4. 

4Q The hybrid loop referred to hereinbefore is obtained 
by removing all terms from the sum of Equation (26) 
exce pt@and@) For example, for low SNR’s 


^- ^ [^VP. sin 

-MV*S[cos 0n, 3 + sin +n ua nis|-f 91 


For high SNR’s the equivalent total phase noise is 
obtained from Equation (31) by replacing the term in- 
volving gg 2 g 3 * by g 3 g3K 2 d n, 3 . In the diffusion approxi- 
mation technique, the coefficient Koo is characterized 

by . 


Mechanization is illustrated in FIG. 5. The operator 
50 Q( ) is simply a multiplication by unity for low SNR’s 
and may be omitted. For high SNR’s, the operator 
would be mechanized as a hard limber, as noted herein- 
before. The subcarrier estimates 5 and — 5 in FIG. 5 
and the low-pass filters can be omitted in the receiver 
55 structure at the expense of additional noise. Such loops 
are of interest in command, low-rate coherent teleme- 
try systems and military applications where the phase- 
error is not constant over the symbol interval. 

Combinations of the data-aided and hybrid loops are 
60 also of interest. When the phase-error is constant dur- 
ing a symbol interval one should take advantage of the 
independence of the noise which is forcing the loop as 
well as the power in the sidebands. In this case the loop 
equation is obtained from Equation (26) by removing 
65 all terms from the sum except the even ones and adding 
The loop equation is then given by 

0= 0 ~ KyF(p)lp [%)+%+•%)+ (36). 
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It is a simple matter to obtain the probability density 
function of the phase-error since one must simply com- 
bine the theory given in the references cited hereinbe- 
fore in connection with the hybrid loop and the data- 
aided loop with the theory given by the inventors in 
“The Performance of Suppressed Carrier Tracking 
Loops in the Presence of Frequency Detuning” and use 
the general theory developed by one of the inventors, 
W. C. Lindsey in “Nonlinear Analysis of Generalized 
Tracking Systems”. A typical mechanization of the 
loop is illustrated in FIG. 6 . As in other cases, the oper- 
ator Q( ) is simply a multiplication by one for low SNR 
and is best mechanized by simply omitting it, and is sgn 
x for high SNR mechanized by a hard limiter. 

Suppressed carrier loops are of interest in practice at 
both the carrier and subcarrier level. Various mechani- 
zations will now be described which render improve- 
ment in such loops when the phase-error is constant 
over the symbol interval. When the carrier is sup- 
pressed, m = P c = 0. For this case, the loop Equation 
(26) reduces to 

$= 6 - KyF(p)lp [%>+^ (37). 


14 

resent the receiver’s estimates of the data subcarrier 
and the data waveforms, respectively; 

„ „\KtK,P 0 sin 2« 

S ® = 0 ‘ { 2 

+VP^[<Kj cos 0»ni — Ki sin 0nu] — nuinul 


where n ul and n lt are respectively the noise processes 
which emerge after separate lowpass filtering of in- 
phase and quadrature phase detections of said signal x 
and the signal %>is the product the lowpass filtered in- 
; phase and quadrature phase detections of said signal x, 
and g, is the d-c gain of said lowpass filtering processes; 
i 

j %= g t K, [ sin<f> + n UI ] 

’where the signal Sgis said inphase phase detection of 
said signal x, neglecting double frequency terms, and 
represents dynamic phase error; 


S®— ffi?a^i a ^ T ^TV'Sf > o-^-^+VP^dT sin 0nuiT 


and the probability density function of the phase-error 
is easily obtained as before. Moments of the means 25 
time to first slip and the average number of slips per 
unit of time can be obtained by using the general theory 
given in the last reference cited. Various mechaniza- 
tions of the loop are possible using the circuit of FIG. 

3 to produce <?(/— T) 

What is claimed is: 

I. In a receiver chan nel fo r a time v aryin g signal, x, 
characterized by x — V 2P C sin <l> + V 2 S X cos <E> + 
n h where =<o „t+6, 0 characterizes modulation due to 

receiver motion or the randomness of said channel, n, 
is a narrowband, “white” Gaussian noise process of 
double-sided bandwidth W ( Hz and single-sided spec- 
tral density N 0 watts/Hz, P c = m 2 P represents power at 
the carrier frequency, 5 = (1— m 2 )P represents the 
power remaining in the modulation sidebands, m de- 
notes the modulation factor, X = 5<^ is a biphase modu- 
lated subcarrier, and 5and d represent the data subcar- 
rier and the data waveforms, respectively, which are 
assumed to be square waveforms, a tracking loop com- : 
prised of a voltage controlled oscillator for generating) 
a time varying reference signal r„= VI/^cos® at an 
output terminal thereof in response to a feedback sig- 
nal at an input terminal, where said feedback signal is 
the sum of one or more feedback signals ^through S© 
at respective points(I)through(2)of said loop excepting 
a signal at poinl®(2)or(§)by itself, or a signal at point 
(g)by itself for low data stream signal-to-noise ratios,, 
and excepting sums of only signals S*jpnd %pr only sig- 
nals S^nd S@for low data stream signal-to-noise ratios, 
said loop including means responsive to said received 
signal and said reference signal for generating said one 
or more feedback signals, where said feedback signals 
are characterized by the following equations, neglect- 
ing double frequency terms: 60 

e„ = /f|[ V^ sin <f>+ VS x cos<f> + «„((,</>)] 

K 2 l Vp7 cos <t> - V 5 x sin<j> n,( (,<£) ] 
where «„ is the output of an inphase phase detector and 
f, is the output of a quadrature phase detector using the 
reference signal r u for inphase phase detection of said 65 
signal x and a 90° phase shifted reference signal charac- 
terized by r, = v 2 K 2 sin for quadrature phase de- 
tection of said signal x and X —Sd is a biphase modu- 
lated subcarrier, 5 and d represent the data subcarrier 
and the data waveforms, respectively, and ^and rf rep- 


+Vs<m a t cos 0n n ,+dT»in o , T ] 


where said estimate 5 of the A data subcarrier is a square- 
wave and it is assumed that 5 = 5 > and the inphase phase 
detected data subcarrier of said signal e„ is filtered by 
30 a lowpass filter of gain g 2 to provide a signal repre- 
sented by g 2 K, [d V""5 cos <f> + n^] which, upon being 
delayed for one data symbol period T and multiplied by 
said data waveform estimate d, delayed one data sym- 
bol period T is multiplied by said signal %>to yield said 
35 signal %>at point®; # * 

%= [ V5 dd sin(£ + d /i„] 

where a represents a data stream estimate for high data 
stream signal-to-noise ratios produced by said bandpass 
filter of gain g 2 in cascade with a bandpass filter of gain 
40 g 3 cascaded with a generator of a function Q{x) ~ sgn 
x implemented as a hard limiter for the inphase phase 
detected and subcarrier detected signal x, and n n is the 
noise process, which is approximately low-pass Gauss- 
ian for the quadrature phase detected and subcarrier 
45 phase detected signal of said input signal x, and said sig- 
nal S 0 is produced by multiplying the output of said 
function generator by the quadrature phase detected 
| and subcarrier phase detected signal of said input signal 
x, and said signal S^is produced for low data stream sig- 
50 nal-to-noise ratios in the same manner, but with said 
generator of a function ,Q(x) = x, where x<<1, in 
which case said signal %>is given by the following equa- 
tion: 

55 SQ=[ffi‘gs 1 KiK 2 (d-/S cos 4>+n aS )(d^/S sin04-n„)] 
=^^K, 2 ir 2 [^|^+<iVS (cos 0 n,» 


-j-ain 0nuj) + Wuanit J 


and said signals Sgand %)at points® ®and@are 
given by 

S®=0j 8 -Ki2G[<2t(<JtVS cos 0 + nuir) < 2 t ( <WS sin 0 +nuT)] 
= gfKt ^ - y 2 * +d r ^ (sin^nuJT 

+cos 0 Wi2t) “b fluJTJiirr^ 
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S@=fcJP[V3(l-2P B (*)) Bin <£-{- c^tWi 2 t] 

^0=t7i^ai^2 2 (VK cos <t>~ ^ii)(VS(1-'2Pe( 0) ) sin £ 

+ <2tWut) 

=^^[VP^S(1-2P e ( 0)) ^+V^rCOB*» UT 
— ^t^uWot — V^( l~ 2 Pe(^)) sin ^nul 


where the signal at point(2)is analogous to the signal at ; 
point(3)with phase error 4> shifted by 90°. 

2. The combination of claim 1 including means re- 
sponsive to said received signal delayed one symbol pe- 
riod T, said reference signal and said reference signal 
shifted 90° for producing one or both of respective de- 
layed inphase and quadrature phase detected and fil- 
tered signals %pnd %>at respective pointsd^nd^given 1 

by 

%>= giK,[ VP, sm<$> + n u ,rl 
•S®= g,X 2 [ y/~F c cos<f> + n tlT ] 
where it is assumed that the phase-error <t> is constant 
over the symbol period T, and the noise process n ulT is, 
orthogonal to the noise processes n ul , n u2 , n u3 , n n , n n 
and n n when the correlation time of the noise is much 
less than T, and said signals at points(Dand®produced 
are added to signals produced as corresponding in- 
phase and quadrature phase detected and filtered sig- 
nals of said undelayed input signal Jt, where the inphase 
phase detected and filtered signal is said signal S^at 
point© 

3. In a receiver chan nel fo r a time varying signal x 
characterised by x = V lP c sin 0 + 2 VS X cos <t> + 

, where X —5 d is a biphase modulated subcarrier, 5 
and d represent the data subcarrier and the data wave- 
forms, respectively, which are assumed to be square 
waveforms, 5 and 3 represent the receiver’s estimates 
of the data subcarrier and the data waveforms, respec- 
tively, and where P c = m 2 P represents power at the car- 
rier frequency, S = ( 1 —m 2 )P represents the power re- 
maining in the modulation sidebands and m denotes the 
modulation, a generic tracking loop, provided to ex- 
1 ploit the principle of nonlinear coherence comprised! 
j. of: ; 

a voltage controlled oscillator fo r g enerating a time! 
varying reference signal r u = Vz cos u>, where) 
<P is the time varying loop estimate of <t> and <t> =; 
0 ) o t + 0, where 0 characterizes modulation due to; 
, recei ver m otion or the randomness of said channel;] 
a summing junction and a smoothing filter coupling 
the junction to a control terminal of the oscillator; 
a 90° phase-shift network for providing a quadrature 
phase reference signal r, = VT K 3 sin <t>; 
two multipliers responsive to the receiver signal and 
the signals r„ and r, for producing quadrature phase 
error signals e„ = xr u and €, = xr,; 
a first low-pass filter of a particular bandwidth and 
gain coupling the signal r u to a point©connected 
to said summing junction; 

a first multiplier having one terminal connected to 
receive the output of said first filter and the output 
of a second low-pass filter of a particular band- 
width and gain to provide a product signal at a 
point(])connected to said summing junction; 
means for demodulating the phase error signal e u by 
a phase estimate of a reference square-wave sub- 
car rier and a t hi rd low-pass filter of a particu lar 


bandwidth and gain for filtering the demodulated 
signal; 

means for delaying this subcarrier demodulated and 
filtered signal a time T equal to a data symbol pe- 
5 riod; 

mqans for multiplying this first delayed signal by 
d( t—T) whered(t) is the time varying estimate of 
the data waveform; 

means for multiplying the output of this last multiply- 
10 ing means by the output of the first filter to pro- 
duce a third feedback signal at point(2)connected 
to the summing junction; 

means for demodulating the phase error signal €, by 
a phase quadrature estimate of the reference 
15 square-wave subcarrier and a fourth low-pass filter 

of a particular bandwidth and gain for filtering the 
phase quadrature demodulated signal; 
means for delaying this subcarrier phase quadrature 
demodulated and filtered signal a time T; 

20 mqans for multiplying this second delayed signal by 
d(t—T) to produce another signal at a point(S)con- 
nected to the summing junction; 
means for multiplying the output of this last multiply- 
ing means by the product of the multiplying means 
25 of the first delayed signal and a(t—T) to produce 
yet another signal at a point ©connected to the 

summing junc tion; 

means for multiplying the output of the second low- 
30 pass filter and the output of the penultimate multi- 
plying means to produce a signal at a point(7)con- 
nected to the summing means; 
fifth and sixth low-pass filters of particular bandwidth 
and gain connected to the outputs of respective 
35 ; third and fourth filters; and 

a multiplier having its output terminal connected to 
a point© one input terminal connected to the out- 
j put of said sixth filter and another input terminal 
connected to the output of said fifth filter by an op- 
40 erator which provides a function approximately 
equal to tanh x, where x is the output of said fifth 
filter; 

I wherein the gain of said filters may be selectively set 
to zero to effectively remove signals at points (T) 
45 through © to provide a desired combination of 

feedback signals to said summmg junction. 

4. The combination of claim 2 adapted to exploit 
sideband power in applications where phase error can 
be assumed to be constant over several data symbol in- 
50 tervals, by providing additional feedback signals at 
points(§)and(9)using a delay means of a delay time T 
coupling said receiver input signal x to two additional 
multipliers, one receiving the reference signals r„ and 
the other receiving the reference r, for inphase and 
55 quadrature phase detection of the delayed input signal, 
and using separate low-pass filters of particular band- 
width and gain coupling the outputs of said additional 
multipliers to said points(S>and® means for connecting 
the filtered output of the inphase error signal thus pro- 
60 duced at point(g)to said signal %>at point© and means 
for connecting the filtered output of the quadrature 
phase error signal thus produced at point(5i wherein 
the gain of said separate low-pass filters may be selec- 
tively set to zero to effectively remove signals at points 

65 (§)and(9)from said tracking loop. 

* * * * * 



